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Magnetic focussing of electrons and holes in the presence of spin-orbit interactions 
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In this work we theoretically investigate transverse magnetic focussing in two dimensional electron 
and hole gasses with strong spin orbit interactions. We present a general result for spin orbit 
interaction with singular winding numbers in the adiabatic limit. We then present results for 
systems with two spin orbit interactions of different winding number appear, using the concrete and 
experimentally relevant case of an applied in-plane magnetic field in hole systems with Rashba type 
interactions. We predict that the application of a large in-plane field is found to have a strong effect 
on the magnetic focussing spectrum. 

PACS numbers: 72.25.Dc, 71.70.Ej, 73.23.Ad 


I. INTRODUCTION 

The dynamics of charge carriers in spin-orbit coupled 
systems is a vital area of investigation for the extremely 
active field of spintronics. Controlling and manipulating 
the flow of electrons and holes serves as the foundation 
of an entire class of spintronic devices, most notably the 
Datta-Das spin transistor and its progeny^. Experimen¬ 
tal studies of spin-orbit systems, with an eye to these 
potential applications, are extensive. One particularly 
fruitful experimental technique is transverse magnetic fo¬ 
cussing (TMF), which involves the mesoscopic transport 
of charge carriers from an injector quantum point con¬ 
tact (QPC) to a collector QPC, in a two dimension (2D) 
charge gas, through which the charge carrier can prop¬ 
agate ballistically on a scale of tens of microns. The 
charge carriers are “focused” by out of plane magnetic 
field B = Bz as illustrated in Fig. [T] The two pos¬ 
sible trajectories shown in Fig. [T] correspond to two 
possible spin polarisation. The figure illustrates that 
the problem under investigation is how the spin influ¬ 
ences the semiclassical long-range orbital dynamics. As 
an experimental technique, TMF has provided impor¬ 
tant insights into the properties of heterostructures and 
Quantum point contacts^. More recently, TMF has been 
applied in Graphene^. Another recent development is 
the use of two dimensional systems with large spin orbit 
(SO) interactions in TMF, the goal to observe spin split¬ 
ting in the magnetic focussing spectrum. Such a split¬ 
ting has been observed in heavy hole systems using zinc 
blend heterostructure^. Simultaneously there has been 
extensive theoretical investigation of TMF with Rashba 
SO interaction which is linear in the particle momen¬ 
tum using a variety of exact numerical^, and strong and 
weak coupling semiclassical approximations^ii. In addi¬ 
tion, polarised photocurrents have been experimentally 
detected using TMF—, and TMF has been used as an ef¬ 
ficient method for detecting density differences in spin 
species in a QPCSiiS. 

While extensive, theoretical investigation so far has 
been limited to the linear in the particle momentum SO 
interactions, like usual Rashba interaction. Though of¬ 
ten dominant in two dimensional electron systems, for 



FIG. 1: In magnetic focussing, electrons (or holes) are fo¬ 
cussed with weak magnetic field from an injector quantum 
point contact (QPC) to a collector QPC, at a distance L. 


heavy hole based systems the spin-orbit Rashba-like in¬ 
teraction is cubic in momenta. Particular crystal lattice 
orientations with respect to the 2D heterstructure in zinc- 
blend semiconductors result in cubic in momenta Dressel- 
haus SO interactioiiii. Application of an in-plane mag¬ 
netic field to a heavy hole zinc-blend semiconductor het¬ 
erostructure results in an effective SO interaction which 
is quadratic in the particle momentum—. Furthermore, 
most 2D systems contain SO interactions of different or¬ 
ders in momentum simultaneously. Such a combination 
can suppress or enhance the spin splitting, depending on 
the growth orientation and confinement shape. Such par¬ 
tial, or complete compensation of spin orbit interactions 
forms the basis of many proposed spintronic devices. 

In light of these considerations, in this work we pursue 
a theoretical investigation of TMF for a wide variety of 
SO interactions which scales as some power of the par¬ 
ticle momentum. Moreover, we consider a case when 
SO interaction with different powers of momentum act 
simultaneously. We assume that the focusing magnetic 
field is weak and hence the evolution of the particle wave 
function is adiabatic. 

The paper is structured as follows. In section m we 
present theory for SO interactions with a single winding 
number, and determining the peak position in a mag¬ 
netic focussing experiment. In section Hill we extend this 
to the case of two SO interactions acting simultaneously, 
and determine the change to the magnetic focussing spec¬ 
trum. Finally, in section HVl we briefly examine the ques- 
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tion of injection via the QPC with a greatly simplified 
model, determining the angular distribution. 


II. SPIN-ORBIT INTERACTIONS WITH A 
GIVEN WINDING NUMBER 

Semiconductor heterostructures allow for a great di¬ 
versity of SO interactions. To render our results as more 
general, we consider the following spin orbit interaction 
for a two dimension electron or hole gas (2DEG), 

'^so = i^p1a+ + h.c. (1) 

Here p is the particle momentum, p± = p^ ± ipy, and 
a± = (Ta, ± iUy, are the Pauli matrices describing the 
effective spin 1/2. For electron systems, these matrices 
represent the electron spin, while for holes which have 
the internal angular momentum J = 3/2, the matrices 
describe the two level heavy hole subsystem Jz = ±3/2, 
with light holes, Jz = ± 1 / 2 , having a signihcantly higher 
energy. 

Derivation of the kinematic form of the interaction 
0 is straightforward. In the case of electrons the only 
possible kinematic structure is the Rashba interactioiJ^, 
'Hso oc (n • [p X S]), where S is electron spin and n is 
the unit vector orthogonal to the plane of 2DEG. This 
results in Eq. with n = 1 and with a real coefficient 
7 = 71- 

For holes in an in-plane magnetic field B the only pos¬ 
sible kinematic structure is Hso oc (J • p)‘^{J ■ B). To 
project this Hamiltonian to the Jz = ±3/2 heavy hole 
subspace one has to replace J/. —>■ cr+, —>■ (t_. All 

other powers of J are projected to zero. Hence Eq. ® 
has n = 2 and a complex coefficient 7 = 726 where 
72 is real and proportional to the magnetic field. The 
phase ip is the angle of the in-plane magnetic field with 
respect to the y-axis, = —i?|| sin(/ 5 . By = i?|| cosp, see 
Fig. [H 

The only possible kinematic structure for the heavy 
holes Rashba-like interaction is Hso oc (J-p)2(n.[px J]). 
Again, projecting it to the Jz = ±3/2 subspace with the 
replacements Jj/ —?► ct+, J/. —?► ct_ we come to Eq. o 
with n = 3 and with a real coefficient 7 = 73 , 

An odd n in Eq. o implies the change of sign under 
inversion, therefore these interactions are due the lack 
of inversion symmetry in the semiconductor heterostruc¬ 
ture. In Rashba-like interactions the inversion asymme¬ 
try is described by the vector n. Alternatively, the n=3 
case can originate from the Dresselhaus SO interaction 
which is due to the lack of inversion symmetry in the bulk 
of a zinc blend semiconductor—. In this case the mag¬ 
nitude and the phase of 7 depends on the orientation of 
the heterostructure with respect to crystal axes. 


The coefficient 7 in Eq. o can be presented as 

7 = 7ne“*‘^ 

ip = 0 if n = 1, 3 

ip ^ 0 if n = 2 

ln=ln-r^ , kF = \/2meF , (2) 

txp 

where ep is the Fermi energy (chemical potential). The 
nonzero phase for n = 2 is determined by the orientation 
of the in-plane magnetic field as it is explained above. 
The dimensionless coefficient 7 „ represents the value of 
the SO interaction a,t p = kp'm units of the Fermi energy. 
For the Rashba-like interaction (lack of inversion asym¬ 
metry in the heterostructure) at n = 1,3 the coefficient 
can as large as | 7 „| ^ 0.1 — 0 . 2 , dependent on the het¬ 
erostructure quantum well. For the n = 2, the interaction 
induced by an in-plane magnetic field, in GaAS the value 
of the coefficient can be a large as I 72 I ~ 0 . 02 /resZoi^, 
depending on the particulars of the heterostructure. The 
Dresselhaus interaction results is a relatively small value 
of I 73 I, in GaAs I 73 I ~ 0.01. It is worth noting that in 
this case the phase can be nonzero. 

We can consider the SO interaction as a momentum 
dependent effective Zeeman magnetic field, H(k). From 

Hso = -Ba- (3) 

B = 7 „/c"(— sin(n0 ± (p), cos{n6 -|- </?), 0) 
k = fc(cos0,sin0,O) , 

where 9 is the axial angle in the /c-space. In the absence of 
a focussing out-of-plane magnetic field electrons (holes) 
move in straight lines, and the spin is parallel (or anti¬ 
parallel) to the SO effective magnetic field, 

(XskIXs) =s(-sin(n 6 »-b(/ 5 ),cos(n 6 l-b<p), 0 ) (4) 

i/k^x) oc e Xs , Xs= { j 

where s = ±1 is the spin polarization, ■i/ks is the par¬ 
ticle wave function, and Xs is the spin wave function. 
Note that the effective spin rotates n-times when k ro¬ 
tates only once in the momentum space. For n = 1 this is 
illustrated in Fig. [5^. This is why we say that the Hamil¬ 
tonian o has a given winding number. A combination 
of Hamiltionians with different values of n obviously vio¬ 
lates this property. It is worth noting that for n = 1 the 
winding number represents the rotation of the real spin. 
However, this is not true for higher values of n. The 
point is that for n = 2, 3 the effective spin a describes 
phases of the heavy hole components Jz = ±3/2. The 
expectation value of (cr) is not equal to the expectation 
value of J. For example for n = 2 (J) does not change 
its direction, it is directed along By, so it does not wind. 
Nevertheless, for our purposes the effective spin matters, 
and therefore the classification by the winding numbers 
makes sense and is important. 
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FIG. 2: The spin orientation for linear in momentum Rashba 
interaction in the k space, (a) In the absence of a focusing 
magnetic field. A particle has a definite momentum, so differ¬ 
ent points in the k-space correspond to different particles, (b) 
A particle in the focusing magnetic field. Both the particle 
momentum and the particle spin evolves in time. Importantly, 
there is a nonzero z-component of spin. 


following equation 

i^U{t)xsO = UU {t)xsO (8) 


= -Bo ■ (TXsO (9) 

Bq ■ cr = ■ crU + 

Bo = (0,7„A:”, -b ^nuj) . 

Solution of Eq. © is straightforward since Bo is time inde¬ 
pendent. There are two eigenstates with energies ±|So| 
corresponding to s = ±1. In a superposition of these 
states the spin oscillates around direction of Bo- In the 
eigenstates the effective spin has the following values 


{<Jx) = 0 


i'^y) 

{(^z) 


luk^ 

'\Bo\ 

{qBbBz -b nuj)/2 

\Bo\ 


( 10 ) 


It is worth noting that the term proportional to nuj in 
(ctz) is a direct consequence of the Corialis force. Eq. 
(flOl) gives the effective spin in the corotatine frame. In 
the laboratory frame this gives 


The full Hamiltonian, in a transverse magnetic field is 



y7r!(cr+ -b h.c. 


(5) 


where tt = p — eA with A being the vector potential 
corresponding to the focusing field B B = V x A. 

Besides the vector potential we include Zeeman coupling 
to the focussing field, g is the g-factor for the electron 
(or hole). Hereafter we include the Zeeman coupling in 
the effective SO magnetic field, 


B 7 „fc”(- sm{ne+ip), cos{n9+(p),0)+l-ggBBz{0, 0,1) , 

( 6 ) 

yielding a finite z component of B. 

When a weak transverse magnetic held is applied, the 
trajectory of the charge carriers is curved, rotating the 
states through some angls 9{t) = ojt, after some time t, 
k —7> k{cosujt, sinojt, 0). For such a curved trajectory, the 
SO held is no longer constant, but rotates in time, see 
Eq.Q. It is conevenient to perform transformation to 
the reference frame corotating with the SO held 


Xs = UxsO (7) 

[/(t) = 6i = wf. 

Substitution of o in Schroedinger equation shows that 
the spin wave function in the corotating frame obeys the 


((T+) = 

(CTz) = s 




1^01 

{gPBBz -b nuj)/2 

\Bo\ 


( 11 ) 


The effective spin evolution in the laboratory frame is 
illustrated in Fig. [5 ]d. 

Now let us look at the spatial dynamics of the particle 
in the laboratory frame. We are interested in a semiclas- 
sical wave packet propagating along a trajectory. To ap¬ 
proach the problem we use the Heisenberg representation, 
where the wavefunctions are time independent, while the 
operators evolve in time. The Heisenberg equations of 
motion with the Hamiltonian ([5]) read 


f+=v+= ilT-L, r+] = — -b nixTr'^ ^a+ 
m 

7f+ = i['H,TT+] = —ieBzV+ . (12) 


The second of these equations represents the usual 
Lorentz force, interestingly it is not inhuenced by the SO 
interaction. On the other hand the first equations which 
gives relation between momentum and velocity depends 
on the SO interaction. In the semiclassical limit we have 
to replace the operators 7 r_|_, u_|_, and in Eas. dT^ by 
their expectation values. The expectation value of spin 
is given by Eq. m- Hence, in the semiclassical limit 
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Eqs. m are transformed to 


V+s = VFse “ 

TT+s = 

kps A 1 - 

VFs = - 1 - 

m \ 2 

iuskps = imWcVFs 
-eB. 


UJc = 


> 0 . 


(13) 


Note that while in Ea. (ll2l) velocity and momentum are 
operators, in Eo. lfOl) these are usual numbers. Note also 
that the Fermi momentum and the Fermi velocity depend 
on spin. When transforming from (IT^ to (IT^ we assume 
that the SO interaction is small compared to the Fermi 
energy, 7 „ <C 1 , we expand in powers (sy^) keeping only 
the leading term. Note that (s 7 n)^ = 7 ^ is independent 
of spin. Therefore Ea. (IT3l) is robust, the neglected spin 
dependent term is of the third order, ^ On the 

other hand we assume that the SO interaction is suffi¬ 
ciently strong to justify adiabatic approximation for spin 
dynamics, 7 „ ^ uiclep- From Ea. dT^ we immediately 
find that the frequency of cyclotron motion depends on 
spin 


LOs = Wc 



(14) 


solutions does not use explicitly the out-of-plane compo¬ 
nent of spin (cTz) presented in Eg. lfTTl) . the Coriolis force 
remains “under the carpet”. In Appendix (TAl) we present 
an alternative solution for n = 1. The alternative solu¬ 
tion is presented in terms of the velocity operator and 
its time derivative. The alternative solution is essentially 
based on (cr^) and the role of the Coriolis force is explicit, 
with the final answer identical to Eqs. (HI and (HU). 


III. SPIN-ORBIT INTERACTIONS WITH 
MIXED WINDING NUMBERS 

In this section we consider the general SO interaction. 
The Hamiltonian is 

n = ^-B-rT. (18) 

Again we assume that the spin dynamics are adiabatic, 
the SO held is everywhere much larger than the cy¬ 
clotron frequency, \B\ ^ lOc- Adiabaticity is the essential 
physical approximation. Besides that we assume that 
\B\ cf- This assumption is not essential, but it allows 
us to perform calculations analytically. Practically the 
inequality is always valid. Because of the spin adiabatic¬ 
ity cr —>• sH/|H| and hence (fTSl) is reduced to the classical 
Hamiltonian 


In a focussing experiment particles are injected at en¬ 
ergy equal to the chemical potential ep independent of 
spin. However momentum of the particle depends on 
spin. Substitution of the spin expectation value (HU in 
the Hamiltonian (H)) gives the following momentum and 
velocity 

kps = kp (^1- 

UFs = — + ^S 7 „[l - . (15) 

m \ 2 J 


H., = ^-s\B\, (19) 

where H is a function of tt. Hence Hamilton Eqs. of 
motion are 

dHci TT+ s dB'^ 
d-K- m \B\ diT- 

TTp = iuicmvp , (20) 

where the z-component in B"^ = B^ + By + Bl can be 
safely neglected. Let us look for solution in the form 


Note that both kps and vps are real, this implies that the 
momentum and the velocity are always parallel. Finally 
the focusing distance, see FigUis 

Ls = vps—= 2rc (I + l-s%] , (16) 

Us \ 2 J 


where 

kp 

rc = - 

muc 


(17) 


is the cyclotron radius. Remarkably, when expressed in 
terms of 7 ^ the focussing length is independent of the 
winding number n. It does depend on the spin polari¬ 
sation, and this leads to the double focusing peak. For 
the case of n = 1 , the solution, dUD, was attained by 
Ziilicke^. 

The solution ( 0 , (HI, (HSU has been obtained in 
terms of momentum and velocity, see Eq. (1121) . This 


7r+ = 

v+ = ne*® , ( 21 ) 

where k is real. Expanding in powers of \B\/ep Eqs. (I19L 
(uni) are transformed to 


k = kp [ 1 + s 


2eF 



6 = ujc 


- JB\ ^ kp ^_,,dB^- 

2e_F 2\B\eF dn- 

4 |B|eF V ^ d^+) 


( 22 ) 


Note that v is generally complex while 6 is of course real. 
Complex V with real k indicates that the momentum and 
the velocity are generally not parallel. 

Integration of Eqs.(H2) is straightforward. To be spe¬ 
cific we apply (HU) to the physically interesting case of 
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heavy holes in asymmetric quantum well and in exter¬ 
nal in-plane magnetic field of few to several Tesla. This 
experimental setup is frequently employed for the polari¬ 
sation of the QPC. The presence of the in-plane magnetic 
field leads to the simultaneous action of n = 3 and n = 2 
SO interactions. 

T-Lso = {’JsTrt + 726 “*“^TT^) a+ + h.c. (23) 
In this case 


02 = ^2^3 ^3 _|_ _|_ g *‘^7r2 ) -I- qlTT^TT^ 

\B\ = eF\j3\b{e) 

b{0) = y/l + 2 ( 72 / 73 ) 008(0 -ip) + ( 72 / 73 )^ (24) 


The condition of the spin adiabaticity implies that b{0) 
does not vanish, \j3\b{9) uidiF- Substitution of (l24ll 
in gives the following 
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m 


6 = ujc 



l + s^6(0) 

3I73I a^y 

® 2 6 ( 6 ») 


3 I 73 I ad) - ic{9) 
® 2 b{9) 


a{9) = 1 -h (5/3)(72/73) cos(6» - (p) + d/'di.d/df 
c{9) = ( 72 / 73 ) sin(6» - ip) . (25) 


The velocity is not necessarily aligned with the momen¬ 
tum, due to the spin orbit interaction. Integration of 
these equations results in the following particle trajec¬ 
tory 


9 = + 9o - f ^^d9' 


m 


X iy = 


kp 


mujc 


{i(e*®° - e*®) 



b{9') -H 3i 


c[e') - 

m. 



(26) 


Hete 9q is the angle which momentum of the particle 
makes with the x-axis at the point of injection, t = x = 
?/ = 0, kx = kcos9o, ky = fcsin^Q. We reiterate again 
that generally the group velocity of the particle is not 
parallel to its momentum. In particular, assume that the 
particle is injected along the x-axis, 9q = 0, ky = 0. This 
does not imply that Vy = 0 , one can easily find from 
Eq.dm that Vy = -572 2 ^ sin ^ 0 . 

In principle Ea. ((26|) solves the problem. However, we 
still need to specify what is the initial value of the in¬ 
jection angle 9q. Intuitively it is zero, 9q = 0. In the 
next section we put this statement on a solid ground by 
considering a model of the QCP injector. Below in this 
Section we set 0o = 0. 

Numerical evaluation of Ea. (l26l) is straightforward. 
Trajectories for 73 = 0.1, s = ±1 and zero in-plane mag- 
netric field, 72 = 0, are shown in panel a of FigjH] by 
dashed blue and red lines. These are simple semicircles 




FIG. 3: Particle trajectories in focucing out-of-plane magnetic 
field. Red lines correspond to the spin polarization s = -1-1 
and blue lines correspond to the spin polarization s = — 1 . 
The value of the “Rashba parameter” is 73 = 0.1. Coordinates 
are given in units of cyclotron radius CZl), the injection point 
is at a; = 1 / = 0. Dashed lines in both panels are identical, 
these are simple semicircles corresponding to zero in-plane 
magnetic field. Solid lines account for a non-zero in-plane 
magnetic field corresponding to 72 = 0.05. Solid lines in the 
panel a are trajectories for the field parallel to the direction of 
injection, p = —-k 12. Solid lines in the panel b are trajectories 
for the field perpendicular to the direction of injection, p = 

0, TT. 


as we discussed in Section II. By solid blue and red lines 
in the same panel a of Figl3]we show s = ±1 trajecto¬ 
ries for 73 = 0.1 and 72 = 0.05. A strong dependence 
of the focusing ditance on the in-plane magnetic field is 
evident. The trajectories depend on direction of the in¬ 
plane magnetic field, Figl3^ corresponds to p = — 7r/2, 
the field is parallel to the direction of injection. In the 
panel b of Figl3] we show trajectories for magnetic field 
perpendicular to the injection direction, (^ = 0, tt. The 
dashed lines are identical to that in Figl3^, 73 = 0.1 and 
72 = 0 (zero magnetic field). The solid lines account for 
the magnetic field corresponding to 72 = 0.05. We see 
that for this orientation of the in-plane magnetic field the 
focusing distance is independent of the field. 

Let us take the focusing shift due to the Rashba inter¬ 
action, Eq. m. 


^Lr = srds (27) 

as a reference value. Then the total focusing shift is 

AL = ALail + SB) , (28) 
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FIG. 4: Dependence of the focusing distance on the in-plane 
magnetic field. 5s is the relative field contribution to the 
focusing distance (relative to Rashba). The plots show 5s 
versus the field orientation angle ip for the folowing values 
of the ratio of the field spin-orbit parameter over the Rashba 
spin-orbit parameter 72/73 = 0.1, 0.25, 0.4, 0.5. At sufficiently 
small 72/73 the angular dependence is approximately sinu¬ 
soidal. 



FIG. 5: Left panel: A model of the injector as a ID quan¬ 
tum channel. Right panel angular distribution of injected 
electrons/holes. 


where 5 b is due to the in-plane magnetic field. Plots 
of 5 b versus the field orientation angle ip for = 

0.1,0.3,0.5 are presented in FigH) One can check an- 
alyticaly that at 72 “C 73 the field correction is 5 b ~ 
—TT^ sin ip. 


IV. INJECTION FROM QPC 


The physics of QPCs is incredibly rich, however, for 
the purposes of this work, we are interested in the peak 
momentum at the QPC exit. For this, it is sufficient to 
consider a greatly oversimplified model, excluding all the 
complexities of QPC physics. We consider the injector as 
a ID channel aligned with the x-direction. The width of 
the channel is w and only one transverse mode is excited 
in the channel, see left panel in FiglS] The orbital wave 


function inside the ID channel is 


Ip oc cosxy , —wj^ < y < w/2 

TT 

K = — 

w 


(29) 


We assume that kx <SC x, hence, since the total energy 
is equal to Fermi energy, we conclude that >c k. kp. The 
Hamiltonian inside the channel reads 


H 


2m 
— ep 


X 

2m 

373 ^ -I- 72 cos (^ ) (Ty -I- 72 sin pax 


(30) 


The spin orbit interaction given by the second line in 
this equation determines the spin polarisation inside the 
channel. 

Following the Huygens principle the fcy-momentum dis¬ 
tribution after emission from the injector is determined 
by the Fourier componenlii^ of (|^ . 


dW 

dky 



cosxy 


2 


(31) 


Converting this to angular distribution we find 

2 


dW 

d9o 


oc cos 00 


cos (1 sin 9o) 


1 — sin^ 0( 


0 


(32) 


The most important conclusion is that the distribution 
(l32]plotted in the right panel of Figgis peaked at 9o = 0. 
This what we used in the previous section. One can say 
that this conclusion trivial, but it is not quite. According 
to Ea. (l25l) this implies that the velocity distributions for 
injected particles with s = ±1 are peaked at nonzero 
angles 


9 


max 


3 _ sin p 


(33) 


The derivation thus far takes no account of adiabatic 
opening of the QPC, reflecting the geometry chosen for 
this analysis. A more realistic geometry includes the adi¬ 
abatic opening of the QPC, from the minimum width of 
the constriction, Wmim to the point at which adiabatic 
transport breaks down, Wmax- At the minimum width 
of the channel, ky = TrjWmin- As the transport through 
the QPC is adiabatic, at the exit, the transverse momen¬ 
tum is ky^x=o = p/Wmax, reducing the angular spread 
of the injected electron/holesJi The magnitude of this 
collimation effect depends on the ratio between Wmin 
and Wmax , and hence the details of the experimental de¬ 
vice. In the above geometry, Wmin = Wmax and hence 
kyWmax 1, leading to a very large angular spread, with 
no diffraction fringes. We note that while flaring can add 
significant collimation to the distribntion, it does not ef¬ 
fect peak, which will still be at fcy = 0, with non-zero 
velocity, as presented in (I55|) . 
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V. SUMMARY 

We consider magnetic focusing of electrons/holes in 
presence of a strong spin orbit interaction. The spin or¬ 
bit interaction is considered in adiabatic approximation, 
spin follows the effective magnetic field. We classify spin 
orbit interactions by winding number n, the number of 
rotations of the effective spin along a close path parti¬ 
cle trajectory. First we consider the case of a singular 
winding number and extend previously known result for 
n = 1 (usual Rashba interaction) to larger values of n, 
especially to experimentally relevant cases n = 2,3. We 
show that in the case of a singular winding number the 
particle trajectory is circular with radius dependent of 
the spin orbit interaction. Second we consider the most 
interesting case of a combination of spin orbit interac¬ 
tions with different winding numbers. Particle dynamics 
in this case are significantly more complex. We derive 
general semiclassical equations of motion linearized in 
spin orbit interaction. Using the developed technique 
we consider in details the case of a combination of n = 3 
(cubic Rashba) and n = 2 (in-plane magnetic field). This 
analysis can be relevant to dynamics of holes in semicon¬ 
ductors and provides a general approach to the problem 
of spin orbit dynamics while ever the adiabatic approx¬ 
imation is valid. We predict that the application of a 
large in-plane magnetic field will have a significant effect 
on the magnetic focussing spectrum, and will be very 
sensitive to the field orientation. 
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Appendix A: Alternative solntion of n = 1 case in 
terms of velocity operator 

Our solution in Section II is obtained in terms of mo¬ 
mentum and velocity, this method is technically the sim¬ 
plest one. An alternative method is to work with ve¬ 
locity without referring to momentum. Intuitively this 
method is more natural in the semicalssical limit. Here 
we present the alternative solution for usual, linear in 


momentum, Rashba interaction. Of course final answers 
of both methods are identical 

We start with usual Rashba spin orbit interaction, 

TT^ 1 

--h 77r (cr X z) - -gusB^^az (Al) 

Zm 2 

which corresponds to the n = 1 for Eq. ©• The velocity 
is 


V = — -I- 7i(T X z , (A2) 

m 

and hence up to terms quadratic in the spin orbit intearc- 
tion the energy is 


mv^ 

2 


■^gfJ-BBzdz ■ 


(A3) 


Heisenberg equation of motion for velocity reads 


V = uJcV X z -b {v X z, (T • z} -b gfj,B(B ■ 2)71 O' (A 4 ) 


We note the presence of a spin dependent term in the 
right-hand-side of the second equation, which has been 
considered as a ’spin force’, in analogy with the usual 
Lorentz, force (the first term)^^. The third term is the 
result of the coupling between the spin-orbit component 
of velocity, and the Zeeman interaction with the out 
of plane field. We stress that the ’spin force’ appears 
only within this technique, it does not appear within 
the momentum-velocity technique considered in the main 
text, Eas. (fT^ . (l20ll . In semiclassical approximation the 
’spin force’ term can be decoupled as {v x z, cr • z} —>■ 
2[v X Due to the ’spin force’ the out-of-plane 

spun projection cr^, which was “under the carpet” in the 
main text solution, here plays an explicit role. Substitu¬ 
tion of {(Jz) from Eq. m gives 

V = Wc (|l - |7i) V X z (A5) 

Hence we get a corrected cyclotron frequency identical 
to Eq. (HI with n = I. We note that the additional 
spin dependent term due to the Zeeman effect does not 
result in any finite contribution to v in the absence of 
the rotation induced correction to {crz{- Due to Ea. (IA3l) 
the particle speed is spin independent, |v| = kp/m and 
we arrive to the same answer as that in Section H with 
n = 1. 
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